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Abstract 

We calculate the QCD radiative corrections to the production of the supersymmetric 
scalar partners of quarks in e + e~ annihilation. We include both the standard gluonic 
corrections and the genuine supersymmetric QCD corrections due to quark-gluino loops, 
and allow for mixing between left- and right-handed scalar quarks which leads to the 
possibility that the two final state particles have different masses. The corrections are 
found to be much larger than the ones affecting the production of spin h particles. 
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1. Introduction 



One of the best motivated extensions of the Standard Model (SM) of the strong, electromagnetic 
and weak interactions is Supersymmetry (SUSY); for reviews see Refs. [[[], [|. Supersymmetric 
theories provide an elegant way to stabilize || the huge hierarchy between the electroweak sym- 
metry breaking scale and the Grand Unification or Planck scales against radiative corrections, 
and its minimal version, the minimal supersymmetric Standard Model (MSSM), allows for a 
consistent unification of the gauge coupling constants, in contrast to the nonsupersymmetric 
SM ||]; moreover it offers a natural solution of the cosmological Dark Matter problem ||. 

The MSSM predicts the existence of scalar partners to all known quarks and leptons. Since 
SUSY is broken, these particles can have masses much larger than the masses of their standard 
partners; however, naturalness arguments suggest that the scale of SUSY breaking, and hence 
the masses of the SUSY particles, should not exceed O (1 TeV). So far, the search for SUSY 
particles at colliders has not been successful and under some assumptions, one can only set 
lower limits of (9(100 GeV) for squarks || and below ~ 50 GeV for sleptons |7j on their 
masses. Higher energy hadron and e + e~ colliders will be required to sweep the entire mass 
range, up to O (1 TeV), for the SUSY particles. 

For a precise prediction of the the production rates of these particles, radiative corrections 
[at least those which are expected to be rather large] should be incorporated. In particular, 
because the strong coupling constant is large, the QCD corrections to squark pair production 
must be included. The QCD corrections to squark pair production at proton colliders have been 
derived recently [§]. At e + e~ colliders, part of the QCD corrections to squark pair production 
[the ones due to virtual gluon exchange and real gluon emission in the case where the squarks 
have equal mass] can be adapted from a result obtained long time ago by Schwinger || in scalar 
QED; see also Ref. [IDf. However, in the MSSM, this result is not complete for two reasons. 



First, it is well known that the SUSY partners of left-handed and right-handed massive quarks 
mix [FT] and the mixing will allow for the two mass eigenstates to be non-degenerate; therefore, 
the Schwinger result must be generalized to allow for the possibility that the final state squarks 
have different masses. Second, in SUSY theories, the gluons have also spin | partners, the 
gluinos, which interact strongly with squarks and quarks; the gluino-quark-squark interaction 
will induce a new type of QCD corrections which has also to be included. 

It is the purpose of this paper to provide the complete result for the QCD corrections to 
squark pair production in e + e~ annihilation, allowing for the possibility of having two squarks 
with different masses in the final state, and including both the standard gluonic corrections 
and the corrections due to quark-gluino loops. 

The paper is organized as follows. In the next section, we will first set the notation and 
exhibit the various couplings as well as the tree-level results which will be relevant to our 
discussion. In section 3 and 4, we will give the analytical results for the standard gluonic 
corrections and for the corrections due to quark-gluino loops respectively; results for the equal 
mass case will also be given explicitely. In section 5, we will discuss the magnitude of these 
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QCD effects. A short conclusion is given in section 6. For completeness, we will list in the 
Appendix the scalar one, two and three-point functions which appear in our results. 



2. Notation and tree— level results 

The production of squark pairs in e + e~ collisions proceeds through s-channel photon and Z 
boson exchanges. The interaction of a neutral gauge boson V = 7, Z with squark current 
eigenstates are described by the following lagrangian 

L qqV = -ieA" £ e q q* d, Z» £ (if - e q s 2 w )q* d, ~q t (2.1) 

i=L,R S WCW i=LtR 

with J| L = ±1/2 and e q the weak isospin and the electric charge [in units of the proton 
charge e] of the squark, respectively; s 2 ^ — 1 — c 2 ^ = sin 2 ^. As previously mentioned, the 
supersymmetric partners of left- and right-handed massive quarks will mix [jllj, the mass 
eigenstates q± and g 2 being related to the current eigenstates qi and q R by 

qi = q L cos 8 q + q R sin 9 q , q 2 = -qi sin 6 q + q R cos 6 q (2.2) 

The mixing angle 9 q is proportional to the mass of the quark q. In the case of the supersymmetric 
partners of the light quarks, mixing between the current eigenstates can therefore be neglected. 
However, mixing between top squarks can be sizeable and allows one of the two mass eigenstates 
to be much lighter than the top quark. Bottom squark mixing can also be significant if the 
ratio of the vacuum expectation values of the two Higgs fields which give separately masses to 
isospin up and isospin down type fermions is very large. After the introduction of this squark 
mixing, the coupling between a gauge boson V and two squarks qi and qj with i,j = 1,2 are 
given by [the directions of the momenta are shown in Fig. la] 

a xl = 2(2/f cos 2 9 q - 2s 2 w e q ) , a 22 = 2(2 if sin 2 6 q - 2s 2 w e q ) , a 12 = a 21 = -2/f sin 29 q (2.3) 
In our convention, the couplings of quarks to the photon and the Z boson are given by 

T^ = -*ee q r , I- z = -— nvq - agl5 ) 

2I 3 q L -4s 2 w e q , a q = 2lf (2.4) 
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When discussing QCD corrections we will also need the squark-squark-gluon and squark- 
squark-gluon-electroweak gauge boson interaction lagrangians which read M 



C gqq = -ig s T a g£[? x dpq 1 + g 2 * ^ g 2 ] 

(2.5) 



c w s w 



i=L,R 



3 



Finally, we need the squark-quark-gluino interaction lagrangian which, in the presence of 
squark mixing, is given by 



C 



ggq 



iV2g s T a q [({?i + a^qi + (#2 + a 2 7 5 )g2 ] g a + h.c. 



(2.6) 



where gs is the strong coupling constant, T a are SU(3)c generators and v i: a, are given by 



-{cos0 q 



sin 6^) = a 2 , a x = -(cos#q + sin# g ) 



-«2 



(2.7) 



In the Born approximation, the total cross section for the production of a pair of squarks 
[possibly with different masses when produced through s-channel Z boson exchange] q~i and q~j 
in e + e~ annihilation, including the finite width of the Z boson, is given by 



a B (e + e- 



QiQj, 



KCL .3/2 



GqV e dij6ij 



+ 



ifll + V e) a 



I6c 2 w s 2 w s - Ml 2563^ (s - M 2 7 ) 2 + T%M 



2.8) 



with s = q 2 = (ki — k 2 ) 2 the center of mass energy of the e + e collider and Ay the usual 
two-body phase-space function, 



A 



ij 



2 2\2 a 2 2 



I, J 



m 2 Is 



(2.9) 



In the case where the mixing between squarks is neglected [as is the case for the supersymmetric 
partners of light quarks and for stop and sbottom quarks if the parameter which describes the 
strength of the non-SUSY trilinear scalar interactions A t) b and the SUSY Higgs(ino) mass 
which also enters trilinear scalar vertices /i, are set to zero] the total cross section for left- and 
right-handed squark pair production simplifies to PUL O 



a B (e + e 



qiqi 



■01 



2 ^Q^G^Qi ^ 

e * ~ 16c 2 w s 2 w s - M 2 



+ 



256s 



(s - M 2 ) 2 + T 2 M 2 



(2.10) 



with 



v q + a q 



■12 



J QR 



v q - a q 



(2.11) 



Note that, as it is expected for the production of scalar particles in e + e~ annihilation, the cross 
section is proportional to the third power of the velocity of the final squarks, (3 q = (1— Am 2 q / s) 1 ! 2 , 
and therefore is strongly suppressed near threshold. The angular distribution is also typical for 
spin-zero particle production, it is given by [6 is the scattering angle] 



da 



B 



dcos^ 



fe+e" 



3 
4 



-sm 2 6a B ( 



e e 



mi, 



(2.12) 



In order to include QCD corrections at first order in as, one needs to consider the diagrams in 
Fig. lb-le, the contributions of which will be discussed in the next two sections. 
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3. Gluonic corrections 



The first set of 0{ag) contributions is due to the standard QCD gluonic corrections. They 
consist first on the interference between the contribution of the tree-level diagram Fig. la and 
the sum of the vertex correction Fig. lb and the squark wave function renormalization Fig. lc, 
where a gluon is exchanged in the squark legs. These two contributions are separately ultraviolet 
divergent and these divergences are regulated using the dimensional regularization scheme. The 
sum of the two contributions is, as it should be, ultraviolet finite but it is infrared divergent; 
this last divergence is regulated by introducing a fictitious mass m g for the gluon. One has also 
to include the contribution of the bremsstrahlung diagrams where a gluon is emitted from the 
external squark lines, Fig. Id. This contribution is also infrared divergent, but the sum of the 
virtual and real corrections is infrared finite as expected. Summing all corrections, the cross 
section at 0(as) can be written as 

3 7T 

where corresponds to the contribution of the virtual exchange of gluons and A^ to the 
contribution of gluon emission from the final state squarks. 

3.1 Virtual Corrections 

Summing the vertex correction and the squark wave function renormalization contributions, 
which as previously mentioned are separately ultraviolet divergent, one obtains for the interfer- 
ence between the tree level diagram and the ones of Fig. lb-lc [after normalizing to the Born 
term and factorizing out 4as/37r] 

,2 ,,2 , \l/2 



a(e + e -> Mj) = cr B (e + e -> Mj) 



Ay = A5 + Ag (3.1) 



\V i r, , 1 /i 2 2m 1 A*i A*j \j 

A ij = l °s—r-^ + TTj2^-^-^)^g- — \ — 

H'g \j 1 — /ij — lij — A i? - 

,2 ,,2 , \V2 



— '1 " A ~ fi) log 1 \ /Xj + V log/, 2 , + Fl I (3.2) 



with the scaled variable ji 2 g = m 2 g / s where m g is a fictitious gluon mass introduced to regularize 
the infrared divergence. The ultraviolet and infrared finite function is given by 

Fl = i[log 2 (l - Vl ) - \og 2 (- yi ) - log 2 (l - y 2 ) + log 2 (-// 2 )] + 2 log log( yi - y 2 ) 

+ log(-yi) log(-y 2 ) - log(l - yi) log(l - y 2 ) + 2Li 2 f — ^— ) - 2Li 2 f J^^l") (3.3) 

with Li 2 the usual Spence function defined as Li 2 (x) = — Jq 1 dtlog(l — xt)^ 1 and the variables 
2/1/2 given by 

2/1/2 = i(l+/*?-/*5±A}/ 2 ) (3.4) 
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Note that the term in the second line of eq. (3.2) [up to a coefficient (1 — /if — fi])/ s] is just the 



Passarino-Veltman [13, 14] scalar three-point function in a particular case 



C (k 1 ,k 2 ,s,m qi ,m g ,m 



1/2 



V 



qj, 



1/2 



2 fog ' \ ' J ' 13 log fi g + F l3 



Ifiifij 



the general form of which is given in the Appendix. 
3.2 Real Corrections 



(3.5) 



In terms of the momenta of the particles in the final state [k is the four-momentum of the 
gluon], the amplitude squared of the process e + e~ — > q^g is given by 



e e 



2„4 



3s 



4 S v 



CqVe&ij&ij 



16c 



+ 



256s 



+ 



s 2 fi 2 



+ 



2s 



+ 



w c w 
2s 



+ 



- Mlf -} 

'ifii+ffi 



r|Mf 



(h-k) 2 (k 2 -k) 2 h-k k 2 -k (fci-fc)(fc : 



6) 



Note that, as a consequence of gauge invariance, the amplitude for the sum of the three diagrams 
of Fig. Id, should vanish when multiplied by the four-momentum of the gluon; this provides a 
good check of the calculation. 

The integrals over the three-body phase space in the general case where the masses of the 
two squarks are unequal [and where both the soft and hard bremsstrahlung contributions are 
added up] can be found for instance in Ref. |L5|. Again, after normalizing to the Born term 



and factorizing out 4as/37r, one obtains for the real corrections Af* with i ^ j 



A 



+ 



+ 



1 - /if - fi] 



A 



1/2 



2 log 2 A - log 2 Ax - log 2 A 2 + 2Li 2 (l - A 2 ) - Li 2 (l - A 2 ) - Li 2 (l - A 2 ) 



,3/2 



-A J/ 2 (l + fi 2 + fi 2 ) + /i 2 log A 2 + fi) log Ai + fi 2 fi 2 log A 



1 + 2/i 2 
+ 4 + — ^logAa 



A 



1/2 

y 



1 + 2/i 2 



A 



1/2 

y 



log Ai + 



(1 - /i 2 - /i 



A 



1/2 



logAo-1 



A 2 - 

log 



+ 



(2 + fi 2 + fi] 



A 



1/2 

y 



logA (3.7) 



with 



An 



2ftifLj 



'i - A - A + x 



v 1/2 ) 
\j ) 



A 



1/2 



2fi j/i 



;i T A ± /i 2 - Af) 



(3.8) 



Note that A^ has the same infrared divergence as A^ but with the opposite sign so that, as it 
should be, the sum of the two contributions = A]^ + A^ is infrared finite. 
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3.3 Equal squark mass case 



In the case where the two scalar particles in the final state have equal masses, the previous 
expressions simplify considerably. In terms of the velocity of one of the final squarks one would 
have for the virtual and real corrections 

1-/3 2 



A v 



1 + /3 2 
P 



1 + /3 2 1-/3 
log 



Li 2 



2/3 
.1 + 0, 



1 + /3 
1 



+ - log 2 : 

4 S 1 + /3 



2 + log 

-/9 



1 + /3 2 1-/3 
log- 



+ log (3 log ■ 



/3 

/3 + 3 



1 + /3 



(3.9) 



A 



log/4, 



+2 log ■ 



1+/? 2 , 1 
1 H tttt- log 



P 



P 



+ log- 



4 

-P 2 



2(3 ° 1 + /3 

- 1 + P 2 

+ 



P 



3Li 2 



3 + 7(3 2 3 - 6/3 2 
_ 1 £_ 

2/3 2 
1-P\ 



1 + 



+ 2Li 2 



4/3 3 

/ 1-/3 N 
, 1+P, 



P\ 1 
— log 



P 



41og/3 
log 2 



4/3 2 



i 

log TTp +log 



2/3 . 1-/9 , 91 4/3 
log ^ + 2 log ■ 



1 + /3 

7T 2 5 , 2 1 - P 

3 I " ° 1 + /3 
1-/3" 



1+/3 °l+p 



log 



1 + /3) 2 °l + /3 



(3.10) 



Adding these two contributions, one obtains the total gluonic QCD corrections to squark pair 
production in the equal mass case 



3 1 + /3 2 , 1- 
2— + 3l ° g — 



P 2 



4Li 2 



1-/3 



2Li 2 



41og/3 + 
1-/3' 



10/3 2 



4/3 3 



5/3 4 1 
log 



1 + P 



P + 1+P 2 



P 



+ 2 log /3 log 



1-/3 



.1 + /3/ V ! + 1 + /3 

This result is in agreement with the one obtained by Schwinger 



3 log log 



(3.11) 



2 ^°l + /3_ 
for scalar QED, once a 



misprint in his eq. (5.4-132) has been corrected; see also Ref. [10 



Let us study the behaviour of this correction in the two limiting situations, /3 — > 1 or 0. In 
the limit where the squark mass is much smaller than the center of mass energy, /3 — > 1, the QCD 
correction approaches the value Ay — * 3, i.e. it is four times larger than the QCD corrections to 
the production of massless spin-1 quarks. In the opposite limit (3 — > 0, i.e. for squark masses 
near the production threshold, one obtains a correction Ay — ► 7r 2 /(2/3) — 2, which exhibits 
the well-known Coulomb singularity near threshold; in this case the perturbative analysis is 
no longer reliable and one has to take into account non-perturbative effects as discussed in 
Ref. |§. 

Finally, let us note that the expression eq. (3.11) can be interpolated by the Schwinger 
formula Pi 



A^^-^(l + /3)(vr 2 -6) 



(3.12) 



which, up to an error of less than 2%, reproduces the exact result given in eq. (3.11) 
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4. Gluino corrections 

Let us now discuss the supersymmetric corrections due to quark-gluino loops. They consist of 
the vertex correction where the partner quark and a gluino are exchanged, Fig. le, and of the 
contribution of the squark wave function renormalization, Fig. If. In the general case where 
squark mixing is present, leading to two final particles with different masses, the expressions of 
these contributions are rather complicated. 

4.1 Vertex Correction 



The contribution of the quark-gluino vertex correction to the coupling of the Z boson to a pair 
of squarks q^ can be written as 



where 



13 4cwSw 



O 7T 



(4.1) 



v q (ViVj + ajOj 



) [<f [(2 



m 2 g 



2mi 



m 2 



<ij 



+ W [{2m 2 g + 2m\ + m\ + m 2 q .)C^ + (m\ - mJ.)C^ + 2m 2 f% + B°(s, m 2 q , m 2 )] } 

+ aqivihj + OiVj) {<f [{2m 2 g - 2m\ - m\. - mJ.)Cj + (m 2 q . - mJ.)C£ 

+ P [(2m; - 2mJ + m\ + mJ)C7j + (m* - m%)C~ + 2m 2 g C% + B°(s, m% w 2 )] } 

+ 2m g m q q^ 2v q {v i v j - a^C^ - a q (vidj - afi^C^ 

+ 2m g m q k»v q (v i v j -a i a j )(2C± + Cl) (4.2) 

where, in terms of the scalar two-point and three point functions B° and C° which are given 
in the Appendix, the functions Cfj = C ± (m 2 .,m 2 ^ s,m q ,m q ,m g ) are given by 



c5 = 



■ ~ 1 
m a 



2sA, : 



B°(ml, m 2 g ,m 2 q ) - B»{m 2 qp mlm 2 q ) + {rh\ - m^C% 
+B°(m 2 m 2 m 2 ) + B°(m 2 m 2 , m 2 ) - 2B°(s, m 2 m 2 ) 



+ (2m 2 



2m 2 g 



m 2 

Hi 



(4.3) 



CZ = 



'2m 2 + 2m 2 - s 



ij 



2s\ 



rh qj )C v 



B°(m^m 2 g ,m 2 ) + B (m 2 q .,m 2 g ,m 2 ) - 2B°(s,m 2 q ,m 2 q ) 
+(2m 2 q - 2m 2 g - m\ - m 2 qj )C^} (4.4) 
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Note that while Cjj is symmetric in the interchange of (& and q~j, is antisymmetric and 
vanishes in the equal mass case. Note also that only the terms proportional to k^ = k% + k 2 are 
ultraviolet infinite [only the two-point function B° is divergent], this has to be expected since 
the one-loop induced terms proportional to = k± — k 2 , which are absent at the tree-level, 
cannot be renormalized. 

The contribution to the coupling of the photon can be straightforwardly derived from the 
previous expressions by setting to zero the terms proportional to a q ; a further simplification 
is obtained by noting that Viv 2 + d\a 2 = and discarding the terms proportional to q^ in the 
equal mass case [since we are interested in the interference with the Born amplitude, these give 
rise to k ■ q 



m qi — m q2 terms which vanish in this case]; one obtains 



4 a, 

3 7T 



6T£ 



(4.5) 



with 



ST?- 



k^iviVj + didj) [(2m 2 g + 2m 2 q + m\. + m 2 q .)C^ + 2m 2 g C% + B°(s, m 2 q , m 2 q ) 
+2m g m q (v i v j - d t h 3 ) [k»(2C± + Cj) + 2q f "C i 



(4.6) 



Note that for % ^ j, the contribution ST^ does not vanish alone [because of the terms in the 
second line of the previous expression: since the two masses m qi and rh qj are not equal, this 



leads to k.q 



™2 



^2 



m„, — m 



7^ 0] and this poses the problem of gauge invariance at the photon 



vertex; however, the left-over piece in the amplitude q^STfj 

q^ST% ~ 2{k.q)C+ + (k.q)C° + 2q 2 Cr j -> B%m 2 q . : m 2 g ,m 2 q ) - B° (m 2 ^ m% m 2 q ) 



(4.7) 



is antisymmetric in the interchange of % and j so that the sum q^ST^ + Qn^ 2 i, which enters in 
the physical process, vanishes and the 7(71(72 is indeed gauge invariant. This feature provides a 
good check of the calculation. 



4.2 Counterterms 



One has then to include the renormalization of the squark wave function. In the general case 
where squark mixing is allowed, this renormalization is a bit more complicated than usual; this 
is because the wave-functions of the two squarks are not decoupled. In addition one has also 
to renormalize the mixing angle q . Taking into account the mixing, the renormalization of 
the squark wave functions and the mixing angle 9 q can be performed by making the following 
substitutions in the Lagrangian eq. (2.1) 

(71 -> (1 + SZ n ) 1/2 <7i + 5Z 12 q 2 , <?2 -> (1 + SZ 22 ) 1/2 q 2 + 5Z 21 q 1 
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Under the substitution 6 q — > 9 q + S9 q , the couplings a^- transform as — > a^- + 5a^-, with 

<5an = 2a 12 56 q , <5a 22 = -2a 12 89 q , <5ai 2 = (a 22 - an)4 (4.9) 

The full counterterms can then be included by shifting the couplings in T^- of eq. (4.1) and 
(4.5) by an amount 



with 



Adij 



Aa u = 5a n + a u 5Z u + 2a 12 SZ 2 i 
Aa 22 = 5a 22 + a 22 5Z 22 + 2a 12 5Z 12 

Aa 12 = 5a 12 + -a 12 (SZ u + 5Z 22 ) + a u 5Z 12 + a 22 5Z 21 



(4.10) 



(4.11) 



One has then to choose a renormalization condition which defines the mixing angle 9 q [and 
hence a± 2 through eq. (4.9)]. We choose this condition in such a way that the total quark- 
gluino QCD correction to the vertex Z^qiq 2 vanishes at zero-momentum transfer q 2 = 0; 5a\ 2 
will be then given by 



5a 



12 



4 a s 1 
372 



s(r 12 + r 21 ; 



g 2 =0 



-ai 2 (5Zu + <5Z 22 ) + an8Zi 2 + a 22 5Z 2 i 



(4.12) 



In the on-shell scheme, where the renormalized squark propagators are such that their poles 
are at p 2 = fh 2 and the residues at the poles are equal to unity, the counterterms SZ^ read 



5Z U = Z' u {m 2 qi ) , 5Z 22 = T! 22 (m\ 2 ) , 
£i2«) s „ E 2 i«) 



5Z 



m% - ml 



, sz 21 = 



ml - ml 



(4.13) 



where Sy(m 2 ) and X-;(m 2 ) = dT,u(p 2 )/dp 2 



are given by 



3tt 



+ a l a j )[A°(m 2 ) + A°(mJ) + (m? + m\ - p 2 )B Q {p 2 , m 2 m 2 )} 



and 



+2m g m q (v i v j - ciid^B (p ,m g ,m q ) 

4«5 



(4.14) 



3tt 



+ aj) [B (p 2 } m 2 m 2 ) + (p 2 - m 2 - m 2 )B' (p 2 , m 2 m 2 )} 



-2m q m g {v i -ai)B (p ,rh ,m 



(4.15) 



with A and B° the scalar one-point and two-point functions which can be found in the 
Appendix. 
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4.3 Case of vanishing mixing angle 



In the case where the mixing angle is set to zero, as is the case for the SUSY partners of the 
light quarks and for stop and sbottom squarks when A t> b = At = 0, the situation simplifies 
considerably. One has just to consider separately left- and right-handed squarks, include the 
vertex correction [with equal squark masses] and the squark wave-function renormalization 
[with the ones of qi and g 2 now decoupled]. The sum of these two contributions is finite and 
no mixing angle renormalization is needed. 

The coupling of the photon and the Z boson to a pair of squarks including the quark-gluino 
vertex correction and the wave-function renormalization, will read at 0(as) 



pM 



pM 

1 ZMi 



4 OCS 1 X 7 



1 + — --A 

3 7T 2 



where in the case of the Z boson, the upper sign is for i 
the corrections factor Af is given by 



-iee q (k 1 + k 2 Y 

±a q )(h + k 2 r [l + -^-A| 

1 =L and the lower sign for % = 2 =R; 



4c w s w 



(4.16) 



A z = 2 



~ 2 i ~ 2 i V q "F a q 2 

m a + m -\ ■ — -m 

y y v„ ± a„ q 



C + + 2m g C + B (s, m q , m, 



D/~2~2 2\ I ~ 1 2 ~2\d'/~2~2 2\ 

-B {m q , m g , m q ) - (m q - m q - m g )B (m q , m g , m q ) 

For the correction to the photon vertex, one just has to replace v q by e q and set a q 
previous equation; one obtains 

A/, = A« = 2(m 2 g + m 2 q + m 2 q )C + + 2rh 2 g C + B (s,m 2 q ,m 2 q ) 
-B (m 2 m 2 m 2 ) - (m 2 - m 2 - m 2 )B' (m 2 m 2 m 2 ) 



(4.17) 
in the 

(4.18) 



In terms of the two and three point scalar functions B and C , the function C + in the equal 
mass case simplifies to 



C 4 



1 



B (m q , m g , m q ) - B (s, m q , m q ) + (m q - m g - m q )C 
Note that the correction to the Z-squarks vertex can be written as 

Af = An + 5u , 5, 



(4.19) 



-m 2 C 4 



1 ± a q /v q 



which exhibits the fact that for the superpartners of the light quarks, where one can set m q 
one would have Af = AZ- = A,-,-. 



(4.20) 
0. 
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5. Numerical results and discussions 



In this section we will discuss the magnitude of these QCD corrections, restricting ourselves to 
the case of vanishing mixing angle [and hence, to the case of degenerate final state squarks] for 
which the numerical analysis is simpler. 

First of all, Fig. 2 shows the magnitude of the QCD correction factor due to virtual 
gluon exchange and real gluon emission, as a function of the velocity of the final squarks, (3 q . 
As previously mentioned, A# rapidly increases from the value A# = 3 for massless squarks (i.e. 
(3 q = 1) to reach an infinite value near the production threshold {(3 q = 0). Close to the latter 
value, perturbation theory is no more reliable and non-perturbative effects have to be taken 
into account. 

In Fig. 3 and Fig. 4, we display the QCD correction factor An due to quark-gluino exchange 
as a function of the squark masses and for two center of mass energy values: = Mz (Fig. 3) 
and yfs = 500 GeV (Fig. 4). While for the SUSY partners of light quarks, only one type of 
correction An (in Figs. 3a and 4a) is present, three types of corrections are needed for the scalar 
partners of the top quark: A] { (in Figs. 3b and 4b) for photon exchange and A^ for Z boson 
exchange and for both the left-handed and right-handed squarks (in Figs. 3c,4c and Fig. 3d,4d 
respectively). 

As representative values of the gluino mass, we have chosen fh g = 5 GeV, 100 GeV and 250 
GeV for LEP energies and m g = 5 GeV, 100 GeV and 500 GeV at y/s = 500 GeV. We have 
allowed for the possibility that the gluinos might be very light, fh g ~ 5 GeV. Indeed, present 
experimental data do not exclude this scenario and a recent analysis |L8| shows that gluinos 
with masses in the 3-5 GeV range are still allowed^. In this case, the bounds on squark masses 
derived from searches for events with a large missing transverse momentum at hadron colliders 
|| might be invalidated^, and thus, the only valid bounds on these masses would come from 
the negative searches of squark pairs in Z decays [7]. These bounds, in fact, will possibly be 
altered by the inclusion of the QCD corrections that we computed here, and this motivates our 
discussion of these contributions for LEP 100 energies. 

At LEP 100, the QCD correction factors due to the contributions of quark-gluino loops are 
displayed in Fig. 3a to 3d for three values of fh g = 5 GeV, 100 GeV and 250 GeV. They are 
practically constant and rather small, |A^| < 1, except for the scalar partners of light quarks 
and for fh g = 5 GeV, in which case A ti varies from —4 for m q ~ 10 GeV to —1.2 near the 
production threshold. The reason is that for large internal particle masses, as is the case at 
this energy for the other two values of m g and for the top quark for which we took a mass 
of m t = 175 GeV ||20|| , the correction is proportional to the inverse of the mass squared; this 
ensure the proper decoupling of the amplitude in this limit. 

^^Note that, even this window would be sufficient to allow for a substantial modification of the running of 
the strong coupling constant |19 which is the original motivation for the renewed interest in this scenario. 

2 This is because of the fact that squarks will dominantly decay into light gluinos, which loose a very large 
fraction of their energy in QCD radiation before they decay [iq] , therefore leading to a soft missing transverse 
momentum spectrum. 
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For the production of squark pairs at a future high-energy e + e~ collider with a center of 
mass energy of t/s = 500 GeV, the QCD correction factors are displayed in Fig. 4a to 4d 
for three values of m g = 5 GeV, 100 GeV and 500 GeV. These factors are larger than in the 
previous case, especially for the scalar partners of top quarks. This is due to the fact that the 
top mass [as well as a gluino mass of 100 GeV] is no more much larger than the beam energy 
and the contributions do not decouple yet. However, the corrections are still relatively small 
except in the light gluino window, and in some cases for gluino masses of the order of 100 GeV. 
Note that the dips in some of the curves correspond to the opening of the q — > q + g channel. 

The QCD corrections to the cross section can be readily obtained by replacing in eq. (2.9), 
the bare couplings e q and v qi with their renormalized values of eqs. (4.16-4.18) [to include the 
corrections originating from quark-gluino loops] and multiplying the bare cross section by a 
factor (1 + 4as/(37r)Ajj) as in eq. (3.1) [to include the standard gluonic corrections]. Writing 
as in eq. (4.20), Af = + 5u, one obtains for the QCD corrected cross section normalized to 
the Born term 



Act 



dee 



a 



B (e+e~ 



1 



4 as 

3 7T 



(5.1) 



with 



da; 



x 



e q v e v qt 



16CjySjy s 



M 2 



+ 



e„ 2 - 



e q v e v qi 



(a 2 e + v 2 )v. 



128S^/C^/ 



m 2 ) 2 



16c 



w°w 



s-M 2 z 



+ 



(a 2 e +v 2 e )v 2 qi 



2564/ct 



is - M 2 ) 2 



(5.2) 



On the Z resonance, one can neglect the photon exchange and the interference between 7 and 
Z exchanges, and the expression of Act simplifies to 



Act 



(Tie e 



a B (e+e' 



qiqi 



^(A„ + A?) 

O 7T 



(5.3) 



In the case of the SUSY partners of light quarks, 8*, 
mass energy, one simply has 



and, independently of the center of 



A 



a(e + e 



CT 



a B (e + e~ 



M) - 1 = |^(A„ + A,,) 

> qiqi) 3 7r 



(5.4) 



The deviation from unity of the QCD corrected cross section normalized to the Born cross 
section as a function of the squark masses and for the previous choices of the gluino masses is 
displayed in Fig. 5 and Fig. 6 for center of mass energies of respectively \fs = Mz and ^/s = 500 
GeV. In these figures, we have used for the strong coupling constant the value acs(M z ) ~ 0.12 as 
determined from various measurements in Z decays [^TJ; at an energy of 500 GeV, the coupling 
will run down to as ~ 0.11. 
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As one might expect from the previous discussion, the largest part of the QCD corrections 
is due to the standard gluonic corrections. The effect of the gluino-quark loops at LEP energies 
is important only in the case of the partners of light quarks when the gluino is very light, 
especially for small squark masses where the correction can even flip sign. At 500 GeV, the 
effect of gluino masses of the order of 100 GeV is also significant and in general, reduces the 
size of the QCD correction. 

For large values of the gluino mass, the quark-gluino contribution decouples and only the 
"standard" gluonic contributions are left. In this case, even for squark masses very small 
compared to the center of mass energy of the collider, the QCD corrections enhance the cross 
section by more than 15%. The correction increases with the squark mass, and already for 
(3 q = 0.5, one has a correction of more than 50%. These corrections are therefore much larger 
than the ones affecting the pair production in e + e~ collisions of spin 1/2 particles with the 
same mass f22|. 



6. Summary 

In this paper, we have calculated the QCD radiative corrections to the production of the 
supersymmetric partners of quarks in e + e~ collisions. We have taken into account the mixing 
between left- and right-handed scalar quarks, and this led us to treat the case where the two 
final state particles have different masses. We included both the standard gluonic corrections 
consisting of virtual gluon exchange and real gluon emission (which, in the limiting case where 
the scalar particles are degenerate in mass, have been calculated long time ago by Schwinger 
for scalar QED) and also the genuine SUSY QCD corrections due to quark-gluino loops (which 
do not appear in scalar QED). For both types of corrections, complete analytical expressions 
were given in the general case. 

We have then discussed the magnitude of these QCD corrections in the case of vanishing 
mixing angle. We have shown that the corrections due to quark-gluino loops are important only 
when the gluinos have relatively small masses, and in this case, particularly for the production 
of the partners of light quarks. For large gluino masses, these corrections decouple as they 
should, and only the contributions due to gluon exchange and real gluon emission are left in 
this limit. The latter corrections are very important since, even for massless squarks where 
they are minimal, they enhance the cross section by more than 15%, four times as much as the 
QCD corrections affecting the production of spin 1/2 massless quarks. 



Acknowledgements: We thank M. Drees and P. Ouellette for discussions. 
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Appendix: Scalar Loop Integrals 

In this Appendix we collect expressions that allow to evaluate the loop functions that appear 
in section 4. The scalar one, two and three point functions, A , B Q and Cq are defined as 



A (m ) 
B (s,m 1 ,m 2 ) 
C (m 1 ,m 2 ,m 3 ) 



(2tt M ) 



n— 4 



(2tt / u)"- 4 



d n k 



k 2 — 171% + ie 



d n k 



(k 2 — m 2 + ie)[{k — g) 2 



(2tt M ) 



n-4 



27T Z 



- m 2 + iej 
d n k 



(A.l) 



[(A; — pi) 2 — m 2 + ie] [{k — p 2 ) 2 — m 2 . + ie](k 2 — m 2 + ie) 



Here n is the space-time dimension and [i the renormalisation scale. After integration over the 
internal momentum k, the function A is given by [7^ is Euler's constant]: 



A (m ) = ml [1 + A ] 



2 u 2 

Ai = Ie + log(4vr) + log ^ 

4 — n mf 



(A.2) 



The function _B and its derivative with respect to s, B' , are given by 

I ™2 ^2 ^2 



B (s,m 1 ,m 2 ) 
B' (s,m 1 ,m 2 ) 



m? — 777,9 m? i + — x_ , 
+ A 2 ) + 2 + log -4 + ^— log - 



1 

~2s 



2 + 



1 1 m i 

-log -4 + 



2 (m 2 — m 2 ) 2 — s(m 2 + m|) 



777o 



x + — X- 



log 



with 



x± = s — m 1 — m 2 ± y s 2 — 2s(m 2 + m 2 ,) + (?n 2 — ml) 2 



(A.3) 



(A.4) 



Note that the x± can be complex; however we have ignored the imaginary parts of Bq and 
B' since they are not relevant for us: to next-to-leading order we are only interested in the 
interference between the (real) tree-level and one-loop amplitudes. 



We also need the three point scalar function Co which can be written in integral form as 



C Q {m x , rh 2 , s, mi, rn 2 , m 3 ) 



dy \ dx ay 2 + bx 2 + cxy + dy + ex + / 



(A.5) 



where 



a = 777, , b = s , c 



7 2 2 

a = m — 777, — 777 



q • 



m 1 — m. 



2 • 



/ = m l 
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(A.6) 



The full analytical expression of Co in terms of Spence functions can be found in Ref. [ I4"| . 
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Figure Captions 

Fig. 1 Feynman diagrams for the QCD corrections to the decay of a gauge boson into squark 
pairs: tree-level diagram (la); vertex (lb) and self-energy (lc) standard gluonic correc- 
tions and real gluon emission (Id) ; vertex (le) and self-energy (If) corrections due to 
quark-gluino loops. 

Fig. 2 The QCD correction factor An due to virtual gluon exchange and real gluon emission as 
a function of the velocity of the final squarks (5. The two squarks are asummed to be 
degerate in mass. 

Fig. 3 The QCD correction factors due to virtual quark-gluino exchange as a function of the 
squark mass and for selected values of the gluino mass at y/s = Mz' A^ for the partners 
of light quarks (3a), A^ for degenerate stop squarks (3b) and Af for the left-handed (3c) 
and right-handed stop squarks (3d). The full line is for m g = 250 GeV, the dashed line 
for m g = 100 GeV and the dot-dashed line for m g = 5 GeV. 

Fig. 4 The QCD correction factors due to virtual quark-gluino exchange as a function of the 
squark mass and for selected values of the gluino mass at y/s = 500 GeV: Ajj for the 
partners of light quarks (4a), A] { for degenerate stop squarks (4b) and Af for the left- 
handed (4c) and right-handed stop squarks (4d). The full line is for m g = 500 GeV, the 
dashed line is for m g = 100 GeV and the dot-dashed line for m g = 5 GeV. 

Fig. 5 The deviation from unity of the fully QCD corrected cross section normalized to the Born 
cross section as a function of the squark masses and for selected values of the gluino mass 
at y/s = M z : for the partners of light quarks (5a), and for the left-handed (5b) and 
right-handed stop squarks (5c). The full line is for m g = 250 GeV, the dashed line for 
m g = 100 GeV and the dot-dashed line for m g = 5 GeV. 

Fig. 6 The deviation from unity of the fully QCD corrected cross section normalized to the Born 
cross section as a function of the squark masses and for selected values of the gluino mass 
at y/s = 500 GeV: for the partners od light quarks (6a), and for the left-handed (6b) and 
right-handed stop squarks (6c). The full line is for m g = 500 GeV, the dashed line is for 
fn g = 100 GeV and the dot-dashed line for m g = 5 GeV. 
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